Astronomy Sz Astrophysics manuscript no. schlattl 


February 2, 2008 


(DOI: will be inserted by hand later) 





Quantum corrections to microscopic diffusion constants 

H. Schlattl and M. Salaris 



co 
o 
o 

(N 
oo 



> 

CO 

in 

CO 
CN 

o 

CO 

o 

0,; 

6 

in 
c3 



S 



Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Egerton Wharf, 
Birkenhead CH41 1LD, United Kingdom 

Received / Accepted 

Abstract. We review the state of the art regarding the computation of the resistance coefficients in conditions 
typical of the stellar plasma, and compare the various results studying their effect on the solar model. We introduce 
and discuss for the first time in an astrophysical context the effect of quantum corrections to the evaluation of 
the resistance coefficients, and provide simple yet accurate fitting formulae for their computation. Although the 
inclusion of quantum corrections only weakly modifies the solar model, their effect is growing with density, and 
thus might be of relevance in case of denser objects like, e.g., white dwarfs. 
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1. Introduction 

Basic physical considerations suggest that, in addition to 
convective mixing - routinely included in stellar evolu- 
tion computations - additional transport processes are 
efficient within the stellar interior; they are driven by 
pressure, temperature and chemical abundance gradients, 
and by the effect of radiative pressure on the individual 
ions. These processes are collectively called 'diffusion pro- 
cesses', and their inclusion in stellar evolution computa- 
tions is necessary in order to sat isfy the helioseismol ogical 
constraint for the solar models (Bahcall et aljfl995j) . 

In general, individual ions are forced to move under 
the influence of pressure as well as temperature gradients, 
which both tend to move the heavier elements toward the 
centre of the star, and of concentration gradients that op- 
pose the above processes. Radiation, which causes negli- 
gible diffusion in the Sun, pushes the ions toward the sur- 
face, whenever the radiative acceleration of an individual 
ion species is larger than the gravitational acceleration. 
The speed of the diffusive flow depends on the collisions 
with the surrounding particles, as they share the acquired 
momentum in a random way. It is the extent of these 'col- 
lision' effects that dictates the timescale of element dif- 
fusion within the stellar structure, once the physical and 
chemical profiles arc specified. The most general treat- 
ment for the element transport in a multi compone n t fluid 
associated with diffusion is provided by iBurgerd ' l|l969l 
B69) equations. In these equations the effect of collisions 
between ions is represented by the so-called resistance co- 
efficients, i.e. the matrices K, z, z' , z", whose precise eval- 



uation is fundamental in order to estimate correctly the 
diffusion timescales for the various elements 1 . 

Recent spectroscopic determinations of Fe and Li 
abundances in Galactic globular clust er t urn-off stars dis- 
cussed in, e.g., iGratton et alJ l)200lh or iBonifacio et alJ 
(2002)), have shown that the present standard treatment 
of diffusion is in disagreement with the observed surface 
abundance of these two elements i n stars of globular clu s- 
ters (for the Li problem see, e.g., iMichaud et aljfl984l or 
IVauclair fe Charbonnell995|h In the light of these results, 
it is very important to conclusively assess how much of this 
discrepancy is due to an incorrect treatment of the diffu- 
sion process, and how much is due to competing rotation- 
ally induced or other non-standard macroscopic mixing 
phenomena, which inhibit the efficiency of diffusion. 

In this paper the state of the art regarding the compu- 
tation of the resistance coefficients in conditions typical of 
the stellar plasma is reviewed and a detailed comparison of 
the results from different authors is performed. The effect 
of different resistance coefficients on solar models is exam- 
ined, too. Moreover, we introduce and discuss for the first 
time in astrophysical computations the effect of quantum 
corrections to the evaluation of the resistance coefficients. 
We also provide simple fitting formulae for accurate calcu- 
lations of resistance coefficients including the appropriate 
quantum corrections. In § 2 we compare existing deter- 
minations of the resistance coefficients, and discuss the 
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1 In the original B69 book the elements of the matrix K are 
called resistance coefficients, whereas z, z , z" , which are re- 
lated to the efficiency of the thermal diffusion, have no specific 
names . In this paper we follow the notation of IPaauette et ail 
(1986), who termed both K and z, z and z" as resistance 
coefficients. 
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differences on solar models; in § 3 we determine the ap- 
propriate quantum corrections to these coefficients, and 
conclusions will follow in § 4. Analytical formulae for the 
computations of updated resistance coefficients, including 
quantum corrections, are given in the appendices 2 . 

2. Classical treatment of element diffusion 

The diffusion of elements in a mult icomponent fluid can 
be tre ated either according to the IChapman fe Cowlind 
<ll970l) or the B69 formalism. We are following the latter 
description, which is equivalent to the so-calle d "second 
approximation" in IChapman fc Cowlind lll97(t . It would 
be desirable to use a higher order approximation, because 
an uncertainty of the or der of 10% is introdu ced by us- 
ing only the second one l)Roussel-Duprelir982l ). indepen- 
dent of the accuracy in the resistanc e coefficients. But this 
can on ly be done in the scheme of IChapman fc Cowlind 
dl970h . which becomes very cumbersome for gases with 
more than two components. 

2.1. Existing calculations of resistance coefficients 

Burgers' equations are obtained assuming the gas parti- 
cles have approximate Maxwellian velocity distributions; 
the temperatures are the same for all particle species; the 
mean thermal velocities are much larger than the diffu- 
sion velocities; magnetic fields are unimportant. Burgers' 
scheme invol ves the resist a nce c oefficients K, z, z' and 
z"; following iMason et all l|l967t MMS67), they can be 
expressed in terms of the so-called reduced collision inte- 
grals ^.' s) * llHirschfelder et alJll954l) 3 . according to the 
following relationships: 
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FORTRAN77 routines to compute the improved 
resistance coefficients (including q uantum correc- 
tions) and their implementation into Schlattl's (2002) 
diffusion-constant routine are publicly accessible under 
http : //www. astro . livjm. ac .uk/~hs. 

3 Note the err o r in the table at page 44 of B69: 
iHirschfelder et alts <1954 8jti0^ j) is equivalent to E^f in B69 



and IChapman fc Cowlinsr s (Il970l) 8^ t } (j). 



Table 1. Quantities in this work, which follows basically 
MMS67, and their relation to quantities used by other 
authors. 
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with = \J k&T '/ (4irn e e 2 ) being the Debye-length, Ay = 
4T 4 * the plasma parameter, fiij = miirij /(mi + rrij) the re- 
duced mass, and Zi 1 mi and rii the charge number, mass 
and particle number density of species i, respectively. The 
collisions between the particles of the stellar plasma deter- 
mine the values of the tiff^ integrals; the physics of the 
collisions is specified by some form of the Coulomb inter- 
action which, as a first approximation, can be described 
by a pure Coulomb potential with a long-range cut-off dis- 
tance, usually assumed to be Ad- According to B69, using 
this truncated pure Coulomb potential the resistance co- 
efficients become 
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where Ce is Euler's constant and the alternative signs in 
Eq. <|7J) represent repulsive (+) and attractive (— ) parti- 
cles, respectively. For high densities and thus small val- 
ues of the plasma parameter, the resistance coefficients 
become negative (for an attractive potential already for 
A < 20). Since in stellar plasmas, in particular in case 
of white dwarfs, A is often much smaller than this value, 
more elaborate computations have been performed. 

More accurate collision integrals have been obtained 
by considering a Debye-Hiickel type of potential 



Vij(r) 
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where r is the particle distance. The results are provided 
either in tabulated form by MMS67 or as fitting f ormu- 
lae bv lMuchmorel dl984l M84).llben fc MacDonaldl l|l985l 
IM85) and lPaauette et ail l|l986l P86). In the latter cases 
the integrals are given as functions of 



$y =ln (ln(l 



A|)) 



(12) 



Different notations for the same physical quantities have 
been used by these groups; in Table ^ their relations to 
the present work are summarized. M84, IM85 and P86 
argue that while Ad is an appropriate screening distance 
at low densities, the Debye sphere loses its significance at 
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Fig. 1. The resistance coefficient K (a), z (b), z' (c) and 
z" (d) as computed by B69 (dotted), MMS67 (solid), M84 
(short-dashed), P86 (dash-dotted) and IM85 (long-dashed 
line). The dark black lines indicate the case of repulsive 
potentials, while the brighter ones show the values for at- 
tractive forces (not computed by M84 and IM85). The 
dots represent the actual values computed by M84. 



Fig. 2. The relative difference of K (a), z (b), z' (c) and 
z" (d) between our analytic fits and the values calculated 
by P86. The solid line represents the case of a repulsive 
potential, the short-dashed one the case of an attractive 
potential; the dotted line indicates the line of zero dif- 
ference. Also shown in (a) is the difference between the 
IM85's and P86's results (long-dashed) and IM85's and 
M84's results (dash-dot-dot-dotted line). 



high densities, in which case a more appropriate screening 
distance is the mean interionic distance; thus they sug- 
gest to use in the actual stellar model computations the 
larger of Ad or the mean interionic distance. Whatever one 
chooses for the actual screening distance, its value has to 
be employed as Ad in Eq. J5J to obtain the appropriate 
value for Ay for determining the collision integrals. We 
also remark that, at least for the entire evolution of the 
Sun up to now, Ad has always been larger than the mean 
interionic distance. 

The results of the different groups are compared in 
Fig.n Noticeably, given that they use all the same physi- 
cal assumptions, they all got very similar values, but they 
all differ significantly from Burgers' results at higher den- 
sities (i.e., lower values of $), where his approximations 
are not adequate anymore. M84's formulae for the case of 
a repulsive potential have been obtained by fitting solely 
the values represented by dots in Fig.^ Although the for- 
mulae work pretty well overall, z" deviates considerably 



at $ = 2 from the mean value of 2.6, and the almost linear 
behaviour of z' for — 3 < $ < 1 is only poorly reproduced. 

The supposedly most accurate calculation of the col- 
lision integrals has been performed in P86, where cubic 
splines for 50 equally spaced intervals in <!> are provided; 
the results agree very well with the tabulation by MMS67. 
To obtain more manageable but still accurate formulae, 
we fitted polynomials of at most 5 th order to their val- 
ues for —2.5 < <f> < 3; this enables to compute diffusion 
of elements up to Fe in stars on the main sequence and 
to follow with sufficient accuracy the early white dwarf 
cooling phase. The numerical values for the polynomial 
coefficients are given in Appendix A. 

In Fig. |21 the relative differences between our fitting 
formulae and the values of P86 are shown. For the case 
of a repulsive potential (e-e or ion-ion collisions) the re- 
sults of P86 for all quantities can be reproduced with an 
accuracy better than 2%; for an attractive potential (e- 
ion collisions) the accuracy is still well within 5%, and for 



4 



H. Schlattl and M. Salaris: Quantum corrections to diffusion constants 



0.005 ^ 
0.004 ~ 




1.2 L_, i , i , , , i , , , i , , , i , , 

0.0 0.2 0.4 0.8 1.0 

r / R e 

Fig. 3. a j The relat i ve di fference between the seismic 
model of Bas u et al and solar models computed 

with different resistance coefficients: SI (solid line), S2 
(dotted), S3 (dash-dot-dot-dotted), S4 (short-dashed), S5 
(long-dashed) and S6 (dash-dotted). The resistance coef- 
ficients used in the models are summarized in Tabic [21 b) 
The run of <E> for H+-H+ (solid) and He 2+ -Hc 2+ collisions 
(dashed line) in solar models. 



a large range of $ values it is certainly better than 3%. 
Taking into account that the results of different groups for 
the same physical assumptions, e.g., between IM85 and 
P86 (long-dashed line in Fig. Et): differ by up to 15%, we 
consider our polynomial fits to be sufficiently accurate. 

In main sequence stars the value of 4> does not de- 
viate considerably from 2 for H and He (see Fig. Et>), 
thus constant values for z, z' and z" are usually assumed. 
Moreover, in stellar plasmas the amount of element diffu- 
sion is basically determined by ion-ion collisions, and as 
the differences for attractive and repulsive potential are 
small when 4> ~ 2, the resistance coefficients for a re- 
pulsive potential are often employed for all c ases. For in- 
stanc e, the widely used diffusion routine by l|Thoul et alJ 
1994, TBL) uses the values of B69 for z, z' and z" (Eq.0- 
I10|) while for K the fitting formula of IM85 has been im- 
plemented. This assumption overestimates z and z', but 
underestimates z" , which results in a somewhat too high 
efficiency of the thermal diffusion (see M84). 

2.2. Effect on solar models 

In order to estimate the influence of different choices of the 
resistance coefficients on the solar structure, we computed 
various solar models utilizing the stellar evolution code 
and clement diffusion routine described in Schlattl (2002) 
and references therein. Element diffusion is treated follow- 
ing the scheme by TBL, considering in addition the effect 
of electron degeneracy and partial ionization of elements; 
the latter is implemented by defining a mean charged ion 
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Fig. 4. a) The surface depiction of various elements in 
solar models with different collision integrals from the pre- 
main sequence until now. The line-styles correspond to 
Fig. |3 On the right hand side the change in the overall 
metallicity is shown, b) Like a), but for the centre. 

per element, instead of computing diffusion for each ion 
separately. This treatment of partial ionization is suffi- 
cient for a large range of stellar masses, including the 
Sun. We neglect radiative levitation, which leads to a tiny 
improvement of the theoretical sound-speed profile com- 
pared to the Sun of at most 0.025% at r « 0.6 R & (see 
iTurcotte et ai1ll99S|) . 

The sound speeds of various solar models computed 
using different choices of the collision integrals are sum- 
marized in Fig. [3^, while the variation of He and metal 
abundances at the surface and at the centre of the Sun are 
shown in Fig. The solid line represents a solar model 
(SI, see Table |2J computed with K of IM85 and constant 
values for the z's of B69 corresponding to the values used 
by TBL. This description has been used, for instance, i n 
solar models of iBahcall et al.l l)l998h and ISchlattl l)2002|) . 
Dropping the assumption of constant z's, and using in- 
stead the functions of M84, the sound-speed difference 
to a seismic model becomes about 4 times higher for 
0.3 < r/i? < 0.65 (dotted line in Fig.EK; S2), caused by 
the reduced thermal diffusion efficiency. The latter leads 
also to a diminished surface depletion of He and metals 
(see Fig.0J. When one employs, in addition, M84's values 
for K instead of the IM85 ones, the difference to model SI 
(dash-dot-dot-dotted line) lowers, because M84 computed 
an about 7% smaller value for K than IM85 (Fig. QJi) for 
the relevant $ range, which results in slightly higher dif- 
fusion velocities. 

In solar conditions, M84's value of 2.62 for z" is how- 
ever too high with respect to the exact result; in fact $ 
is about 2.2 in solar radiative regions (see Fig. Et>) , and 
z" has to be between 2.3 and 2.4 (Fig. [Tjl). Applying the 
correct z"($) would further reduce the sound-speed dif- 
ference compared to models computed with the value of 
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Table 2. Source for resistance coefficients used in the solar 
models discussed in the text and typical helioseismologi- 
cal quantities of the models. The column denoted with z 
represents here z, z' and z"; column "r/a" shows mod- 
els where repulsive and attractive potentials are distin- 
guished, and "qu." where quantum effects are included. 
Note, that in case of P86 our fitting formulae for the col- 
lision integrals have been used. The two columns on the 
r.h.s. denote the resulted depth of the convective envelope 
(Rb.c./ Rq) and the surface helium abundance (Y s ) in the 
respective solar models. 
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r/a 


qu. 


Rb.c./R<3 


Y s 


SI 


IM85 


B69 






0.7125 


0.2448 


S2 


IM85 


M84 
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0.2476 


S3 


M84 


M84 






0.7144 


0.2467 
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P86 


P86 
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0.2457 
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P86 
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M84. Indeed, using our formulae for P86's collision inte- 
grals, which gives similar values of K as M84 (Fig. [20, but 
includes the variation of z" , lowers the sound-speed differ- 
ence (compare S3 and S4 in Fig. |2K)- A small additional 
increase of the diffusion velocities is obtained by consid- 
ering different collision integrals for e-ion (attractive) and 
ion-ion or e-e (repulsive) interactions (see also Fig. |3J|, 
which leads to model S5 (long-dashed line in Fig. Ek) ■ 

Typical additional quantities inferred by helioseismo- 
logical methods are the depth of the convective zone 
and the surface heliu m abundance, which a re determined 
to be 0.713 ± 0.001 llBasu fc Antial Il997|) and 0.246 ± 
0.002 l|Basu fc Antialll995l) . respectively. In all models of 
Table apart from S2 these quantities are well within 
the observational limits. Nonetheless, a general trend fol- 
lowing the overall sound speed behaviour can be ob- 
served: The larger the sound-speed difference of the re- 
spective model the shallower the convective envelope and 
the higher the surface helium content. 



3. Quantum corrections 

All computations of the resistance coefficient discussed in 
the previous section have assumed that the collisions are 
dominated by the classical interaction between two point- 
charge particles. For long-range potentials like the Debye- 
Huckel one, quantum effects become important at high 
energies, whereas t he behaviour is classical at low energies. 
Halm ct a^ l|l97lL H71) have computed corrections to the 
classical collision integrals, such that 



lass ) 

ft^ 9) *(T*A* 2 ) ± 



1 



2s +1 



where f2 



:la.ss 



is the classical value of 32 Slqu 9 '' contains 



0.010 F 



0.005 



0.000 

-0.005 

-0.010 
0.2 

-0.0 



-0.4 - 



-0.6 



:(a) 

- 


, , , , | 




















- 






\ 

1 1 






















\ 






X 




r , . , 1 


\ 




10° 


10 6 10 7 



T (K) 

Fig. 5. a) The relative change of K (long-dashed line), 
z (solid), z' (short-dashed) and z" (dash-dotted) for H + - 
H + collisions in the Sun in the sense (quantum - classi- 
cal/classical. For the dark line the quantum corrections 
tabulated by H71 have been utilized, while the grey lines 
were obtained with our fitting formulae (see Appendix B). 
b) As a), but for H + -e collisions. 



change contribution, which has to be dropped for distin- 
guishable particles (the indices ij denoting different inter- 
acting particles have been omitted for the sake of clarity.) 
The upper sign in the term in parenthesis of Eq. I|13l) refers 
to particles obeying Fermi-Dirac statistic, the lower one 
for particles following Bose-Einstein statistics. The spin 
of the particle is denoted, as usual, by s, and A* is the 
dc Boer parameter given by 



A*. = 

1 y 



2/j,ijZiZje 2 \i 



the quantum mechanical correction, and £1 



the ex- 



Tabulated values for T* 2 fife'&* 

are given in H71, 
which we again fitted by polynomials of 5 th order (see 
Appendix A). 

The excellent accuracy of our fitting formulae can be 
checked in Fig. CD where the effect of the quantum cor- 
rections on the resistance coefficient and the thermal dif- 
fusion coefficients in the Sun is shown. The corrections 
increase toward the centre, and especially for the H + -e col- 
lisions they considerably alter K. However, the diffusion 
velocity of H is basically determined by H + -H + interac- 
tions and in that case the corrections for all quantities 
are below 1%. The expected weak influence on the solar 
structure is demonstrated by model S6, where our fitting 
formulae to the quantum corrections of H71 have been 
used. Overall, the diffusion velocities for all elements have 
been increased, more prominently in the centre, where the 
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difference to models S5 (without quantum corrections) is 
the biggest (Fig. |5j). This caused a further reduction of 
the sound-speed difference of model S6 compared to S5 
(Fig. Si). 

It should be mentioned that the spin-dependent term 
in Eq. (|13f) would demand to know exactly the number 
of atoms in each ionization and in each excited state. 
However, the diffusion velocity of all elements is basi- 
cally determined by their collision frequency with the most 
abundant element (H), and thus only for H + -H + colli- 
sion this spin-dependent term has to be accounted for 
(the spin-dependent term appears only in case of indistin- 
guishable particles). Moreover, when quantum corrections 
become important inside stars most of the elements are 
already fully ionized. Thus only very tiny errors are intro- 
duced when, as in this work, only one mean-charged ion 
per element is considered. 

It is interesting to notice that model S6, which includes 
the most accurate collision integrals of P86 accounts for 
the difference between repulsive and attractive interac- 
tions, and includes quantum effects, has almost the same 
sound-speed profile as model SI (cf. Fig.|3t), which con- 
tains the most crude treatment of diffusion among all mod- 
els considered in this work. Also the depth of the convec- 
tive envelope and the surface helium content are nearly 
identical (Table 0). Therefore, just by chance, a very ac- 
curate solar model with respect to diffusion is obtained 
when using the resistance coefficients chosen by TBL. 



4. Conclusions 

In this paper, we have compared resistance coefficients 
computed by various groups for the calculation of atomic 
diffusion constants, and we have discussed the inclusion of 
quantum effects on the otherwise classical determination 
of these quantities. We provide simple analytical formulae 
to implement easily the presently most accurate determi- 
nation of resistance coefficients, i.e., the classical result by 
P86 plus the effect of quantum corrections determined by 
H71. (FORTRAN routines for computing updated resis- 
tance coefficient and the consequent diffusion constants 
are publicly available. 2 ) 

The different classical computations based on a a 
Debye-Hiickel type of potential produce solar sound-speed 
profiles with significant differences when compared to 
the current accuracy of hclioscismological determinations. 
Just by chance, our accurate treatment including quantum 
corrections produces a solar sound-speed profile compara- 
ble to the one obtained using the less accurate z, z' and 
z" coefficients selected by TBL. In this context, we would 
like to add that with none of the descriptions for the col- 
lision integrals - neither with the inclusion of radiative 
levitation - the bump in the sound-speed difference at 
r w 0.65 Rq disappears. So, an additional mixing process 
beyond the formal boundary of the convective zones is still 
a probable candidat e to resolve this discrepancy (see, e.g., 
Richard et a 1 .ifTiiiit?) . 



The introduction of quantum corrections increases the 
efficiency of diffusion with respect to the classical case, and 
their effect is more pronounced for higher densities. It will 
therefore be important to test the effect of our accurate 
resistance coefficients on models of objects denser than the 
Sun, like white dwarfs, and main sequence turn-off stars of 
galactic globular clusters, where also radi ative levitation 
is altering the surface abundance patterns ijRichard et alJ 

Although the accuracy of the diffusion constants could 
be improved considerably when using the correct func- 
tions for the resistance coefficients, there still remains the 
limitation that Burge r s' form alism is equivalent only to 
IChapman fc Cowlin second approximation. This 

causes an intrinsic uncert ainty when using Burg ers' equa- 
tion of the order of 10% I Roussel-Duprelll982ll . which is 



difficult to reduce further. 
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Appendix A: 

In this appendix we provide analytic formulae for 
T* 2 Q( s '*) and those combinations of the indices s and 
t which are needed to compute the resistance coefficients 
K, z, z' and z" using Eq. The polynomial func- 

tions for the classical case have been obtained by a least 
square-difference fitting of the results of P86. 

The best fits for the classical collision integrals could 
be obtained with 



5 

k=0 



for attractive and 

T *2Q( s ,i)* 
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fc=0 



<£(*) fat) 



with 



$ st = ln(ln(l + A 6 ")) 



(A.l) 



(A.2) 
(A.3) 



for repulsive potentials. We dropped the indices i,j for the 
sake of clarity. The numerical values for the coefficients c A 
and c R , and the exponents b are summarized in Table lA~Tl 
Notice that it is sufficient to determine T* 2 f2( s,t ) , as in 
the definitions of the resistance coefficients the factor T* 2 
always cancels out. 

The quantum-mechanical ( "qu" ) and exchange contri- 
bution ( "ex" ) can both be computed according to the fol- 
lowing relationship 



T* 2 n^li*(T*A* 2 ) = _!^/W ex (T*A 



N- 1 
p 



1 



2(1+ P) ' 



with (A.4) 
(A.5) 



In these formulae the quantum and exchange parts of 
the collision integrals are expressed as functions of the 
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Table A.l. Coefficients for Eqs. (fP|l - lfO)) . 



(s,t) 





(1,1) 


(1,2) 


(1,3) 


(2,2) 


c A (0) 


-1.577 


-2.062 


-2.472 


-1.776 




0.6285 


0.5066 


0.4452 


0.8555 


c A (2) 


0.08141 


0.05224 


0.04911 


0.07976 




03769 


03302 


02851 


01 1 98 


C A (4) 


-0.002702 


-0.0005197 


1.605 xlO -5 


-0.002642 


c A (5) 


-0.001587 


-0.001291 


-0.001014 


-0.0004393 


c R (0) 


-1.862 


-2.465 


-2.857 


-1.702 


c*(l) 


2.313 


1.667 


1.386 


1.916 


c R (2) 


-0.1550 


-0.07154 


-0.05136 


-0.08114 


C R (3) 


-0.07188 


-0.03209 


-0.01620 


-0.04868 


& 


0.6339 


0.8228 


0.9231 


0.7807 







q 






1 


2 


3 


dcx(0) 




-0.7074 


-0.7103 


dex(l) 




-1.985 


-2.744 


dex(2) 




0.5233 


0.1870 


dex(3) 




-0.3270 


0.03388 


dex(4) 




0.03920 


-0.06502 


dqu(0) 











dqu(l) 


0.05333 


0.03051 


0.3133 


dqu(2) 


0.01889 


0.01021 


0.03982 


dqu(3) 


0.007114 


0.009180 


-0.003032 


dqu(4) 


-0.0006953 


-0.0007776 


1.388xl0" 5 


3ex 




6.121 


2.285 


3qu 


0.3 


1.0 


0.15 



quantity T*A* . The best fits to the values tabulated in 
H71 have been obtained with 

/-''•/>, ... - yva 

t=0 



I ( c q) (x) = 52$(i)$\x), where 



£ {q) (x) = I ln ( 1 + -9c )><a; ) for C="qu", 
C [ ' \ ln(l + g {q) / x) for C = "ex" 



(A.6) 
(A.7) 



with x — T*A* and the additional constraint that 
$(x) = for x < 0.1. The index "C" denotes the appro- 
priate expression for either the quantum or the exchange 
contribution to the collision integrals. The respective coef- 
ficients are provided in Table lA~2l Note that Ie* = and 
that Eqs. (|A.6|I and HA.7|> are only valid for x < 1 000, 
as no values were provided by H71 for larger x. We sug- 
gest not to extrapolate T*A* 2 beyond this upper bound- 
ary, but to use the values at T*A* 2 = 1 000, which are 

1.9581 



l£> = 0.4913, 7e ( x } 



and I, 



(3) 

qu 



0.4960, 4i } 



1.5053, 1$ 



2.2010. 



Appendix B: 

Here we describe the modifications to be applied to the 
routine by TBL, in order to use any preferred formulation 
of the resistance coefficients. 



The first change involves Eqs. (9)-(10) of TBL, where 
the IM85 definition of K is implemented. They have to be 
replaced with the chosen representation of K . 

The matrix of elements Y st , introduced at page 830 of 
TBL, has to be changed to 



Y st = 3y s t + z'st x s t m t /m a 



(B.l) 



Finally, the following elements of the matrix Ay as 
defined in Eqs. (33) and (34) of TBL become 



Table A.2. Coefficients for Eqs. (|A~fi)l - l|A~7) 



Z)fc# z ikhkXik for j = i + S, 



for j =5 + 1,..., 25 
and j ^ i + S, 



(B.2) 



for i = 1, . 



, S and 



' 2 - 5 Z)fe#j ZpkkpkXpk for j =p, 



2.5 k"pj "^"pj 



for j = l,...,S and j ^ p, 
Y,k 7 ti k pky P k(0.8Zp k x pk + Y pk ) - OAz'ppk 



for j = i, 
0.8z'; q ) 

for j = S+ 1, 
and j ^ i, 



(B.3) 



,2S 



for i — S + 1, . . . , 25*, where p = i — S and q — j — S. Here, 
the same symbols have been used as in TBL, to which the 
reader is referred for their definitions. 
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